ENUMERATION OF INTEGRAL TETRAHEDRA 
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Abstract. We determine the numbers of integral tetrahedra with diameter d up to isomorphism 
for all d < 1000 via computer enumeration. Therefore we give an algorithm that enumerates the 
integral tetrahedra with diameter at most d in 0(d 5 ) time and an algorithm that can check the 
canonicity of a given integral tetrahedron with at most 6 integer comparisons. For the number of 
isomorphism classes of integral 4x4 matrices with diameter d fulfilling the triangle inequalities 
we derive an exact formula. 



1. Introduction 

Geometrical objects with integral side lengths have fascinated mathematicians for ages. A very 
simple geometric object is an m-dimensional simplex. Recently an intriguing bijection between 
m-dimensional simplices with edge lengths in {1,2} and the partitions of m + 1 was discovered 0. 
So far, for m-dimensional simplices with edge lengths in {1,2,3} no formula is known and exact 
numbers are obtained only up to m = 13 [9]. Let us more generally denote by a(m,d) the number 
of non-isomorphic m-dimensional simplices with edge lengths in { 1 , . . . , d} where at least one edge 
has length d. We also call d the diameter of the simplex. The known results, see i.e. [5], are, besides 
some exact numbers, 
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where p(m + 1) denotes the number of partitions [A000041] of m + 1. The aim of this article is the 
determination of the number of non-isomorphic integral tetrahedra a(3, d). 

Besides an intrinsic interest in integral simplices their study is useful in field of integral point 
sets. These are sets of n points in the m-dimensional Euclidean space E m with pairwise integral 
distances. Applications for this combinatorial structure involving geometry and number theory are 
imaginable in radio astronomy (wave lengths), chemistry (molecules), physics (energy quantums), 
robotics, architecture, and other fields, see [3] for an overview. We define the largest occurring 
distance of an integral point set V as its diameter. From the combinatorial point of view there is a 
natural interest in the determination of the minimum possible diameter <i(m, n) for given parameters 
m and n [U 01 [5J [H [TOl EJ 113 EH EH]. In most cases exact values of d(m, n) are obtained by an 
exhaustive enumeration of integral point sets with diameter d < d(m,n). A necessary first step for 
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the enumeration of m-dimensional integral point sets is the enumeration of m-dimcnsional integral 
simplices. Hence there is a need for an efficient enumeration algorithm. 

Another application of integral tetrahedra concerns geometric probabilities. Suppose you are 
given a symmetric 3x3 matrix A 2 with entries being equi-distributed in [0, 1] and zeros on the 
main diagonal. The probability Vn that A 2 is the distance matrix of a triangle in the Euclidean 
metric can be easily determined to be V2 = \- As a generalization we ask for the probability V m 
of a similar defined (m + 1) x (m + 1) matrix A m being the distance matrix of an m-dimensional 
simplex in the Euclidean metric. To analyze the question for m — 3 we consider a discretization 
and obtain V 3 = lim 4 '"jq . 

d^oo a 

Our main results are the determination of a(3, d) for d < 1000, 



Theorem 1. The number a<(d, 3) of symmetric 4x4 matrices with entries in {1, . . . , d} fulfilling 
the triangle inequalities is given by 

. / 1^+42^+1628^ for d = mod 2; 

, ) \ 17^+42 5 ^ + 176 3 rf 2 +675 ford = 1 mod 2 . 

If we additionally request a diameter of exactly d we have 

j 34d 5 -85^+680rf 3 -962rf 2 + 1776rf-960 ford = Q mod 2 ; 
a( ' 34^-85^+680^-908^ + 1722^-483 /or d = ! mod 2 , 

Theorem 2. 

0.090 < P 3 < 0.111, 

and the enumeration algorithms of Section [4] and Section [5] which allows us to enumerate integral 
tetrahedra with diameter at most d in time 0(d 5 ) and to check a 4 x 4-matrix for canonicity using 
at most 6 integer comparisons. 

2. Number of integral tetrahedra 

Because a symmetric 4 x 4-matrix with zeros on the diagonal has six independent non-zero values 
there are c? 6 labeled integral such matrices with diameter at most d. To obtain the number a<(d, 3) 
of unlabeled matrices we need to apply the following well known Lemma: 

Lemma 1. (Cauchy-Frobenius, weighted form) 

Given a group action of a finite group G on a set S and a map w : S — ► R from S into a 

commutative ring R containing Q as a subring. If w is constant on the orbits of G on S , then we 
have, for any transversal T of the orbits: 

teT 1 1 geGseSg 

where S g denotes the elements of S being fixed by g, i.e. 

S g = {s e S\s = gs} . 

For a proof, notation and some background we refer to [6]. Applying the lemma yields: 
Lemma 2. 

/ . d 6 + 9d 4 + 14rf 2 
a<(d,S) = — 
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and 



a(d, 3) = a<(d, 3) — a<(d — 1, 3) 



Qd 5 - 15d 4 + 56d 3 - 69d 2 + 70d - 24 
24 ' 



As geometry is involved in our problem we have to take into account some properties of Euclidean 
spaces. In the Euclidean plane E 2 the possible occurring triples of distances of triangles are com- 
pletely characterized by the triangle inequalities. In general there is a set of inequalities using the 
so called Cayley-Menger determinant to characterize whether a given symmetric (m + 1) x (m+ 1) 
matrix with zeros on the diagonal is a distance matrix of an m-dimcnsional simplex |13j . For a 
tetrahedron with distances Sij, < i < j < 4, the inequality 
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has to be fulfilled besides the triangle inequalities. 

In a first step we exclusively consider the triangle inequalities for m = 3 and count the number 
a<(d, 3) of non-isomorphic symmetric 4x4 matrices with entries in {1, . . . , d} fulfilling the triangle 
inequalities. 
Proof of Theorem^ 

Counting labeled symmetric 4x4 matrices with entries in {1, . . . , d} fulfilling the triangle inequal- 
ities is equivalent to determining integral points in a six-dimensional polyhedron. Prescribing the 
complete automorphism group results in some further equalities and an application of the inclusion- 
exclusion principle. Thus, after a lengthy but rather easy computation we can apply Lemma [l] and 
obtain 



24a<(d, 3) = 3 • 
34<i 6 + 55d 4 + 
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4d 4 + 5rf 2 ' 
12 

136d 2 l „ 
- 6 
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37d 4 - 18d 3 + 20d 2 - 2ld + (36d 2 + 42) [|] 



96 



which can be modified to the stated formulas. □ 
In addition to this proof we have verified the stated formula for d < 500 via a computer enumer- 



ation. We remark that 



IflS and HI tend to ^ = 0.1416 if d 
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to obtain an exact formula for a(d, 3) because the Cayley-Menger determinant 



CMDo 



for dimension m = 2 can be written as 

CMD 2 = -(S ,i + 5 ,2 + <M(<5o,i 
Thus CMD2 < is equivalent to the well known linear triangle inequalities 5o,i + <5o.2 > $x,2j 
<5o.i + ^1,2 > ^0,2 and <5o,2 + £1,2 > ^o.i- Unfortunately for m > 3 the Cayley-Menger determinant 
is irreducible [I] and one cannot simplify (— \) m+1 CMD m > into a set of inequalities of lower 
degree. So we are unable to apply the same method to derive an analytic formula for a(d, 3). 



+ <5q, 2 - £1,2) (£0,1 - S ,2 + 5l,2)(— #0,1 + ^0,2 + Si 2)- 
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8854161 
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639349793 
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4 


27 


369247 


52 


9756921 


77 


69497725 


140 


1382200653 


3 


16 


28 


442695 


53 


10732329 


78 


74130849 


160 


2695280888 
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45 


29 


527417 


54 


11783530 


79 


79008179 


180 


4857645442 


5 


116 


30 


624483 


55 


12916059 


80 


84138170 


200 


8227353208 


6 


254 


31 


735777 


56 


14133630 


81 


89532591 


220 


13251404399 


7 


516 


32 


861885 


57 


15442004 


82 


95198909 


240 


20475584436 


8 


956 


33 


1005214 


58 


16845331 


83 


101149823 


260 


30554402290 


9 


1669 


34 


1166797 


59 


18349153 


84 


107392867 


280 


44260846692 


10 


2760 


35 


1348609 


60 


19957007 


85 


113942655 


300 


62496428392 


11 


4379 


36 


1552398 


61 


21678067 


86 


120807154 


320 


86300970558 


12 


6676 


37 


1780198 


62 


23514174 


87 


127997826 


340 


116862463817 


13 


9888 


38 


2033970 


63 


25473207 


88 


135527578 


360 


155526991341 


14 


14219 


39 


2315942 


64 
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89 


143409248 


380 
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15 
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40 


2628138 


65 


29783292 


90 


151649489 


400 


263399396125 


16 


27421 


41 


2973433 


66 


32145746 


91 


160268457 


420 


336178761892 


17 


37062 


42 


3353922 


67 


34657375 


92 


169272471 


440 


424224122232 


18 


49143 


43 


3773027 


68 


37322859 


93 


178678811 


460 


529820175414 


19 


64272 


44 


4232254 


69 


40149983 


94 


188496776 


480 


655468974700 


20 


82888 


45 


4735254 


70 


43145566 


95 


198743717 


500 


803900006590 


21 


105629 


46 


5285404 


71 


46318399 


96 


209427375 


520 


978079728301 


22 


133132 


47 


5885587 


72 


49673679 


97 


220570260 


540 


1181221582297 


23 


166090 


48 


6538543 


73 


53222896 


98 


232180129 


560 


1416796092768 


24 


205223 


49 


7249029 


74 


56969822 


99 


244275592 


580 


1688540496999 


25 


251624 


50 


8019420 


75 


60926247 


100 


256866619 


600 


2000468396580 



Table 1. Number a(d, 3) of integral tetrahedra with diameter d - part 1. 



d a(d, 3) 


d a(d, 3) 


d a(d, 3) 


620 2356880503873 
640 2762373382787 
660 3221850132593 
680 3740530243895 
700 4323958989350 
720 4978017317882 
740 5708932993276 


760 6523288334629 
780 7428031732465 
800 8430487428682 
820 9538364312059 
840 10759766492473 
860 12103204603044 
880 13577602128303 


900 15192308794063 
920 16957109053082 
940 18882231158104 
960 20978358597822 
980 23256639532080 
1000 25728695195597 



Table 2. Number a(d, 3) of integral tetrahedra with diameter d - part 2. 



Lemma 3. We have a(3,d) G n(d 5 ), a(3,d) G 0(d 5 ), a<(3,d) G n(d 6 ), anda<(3,d) G 0(d e ). 

Proof. The upper bounds are trivial since they also hold for symmetric matrices with integer values 
at most d and zeros on the diagonal. For the lower bounds we consider six-tuples #o,i G [d, d(l — e)), 
5 ,2 G [d(l - e),d(l - 2s)), 5 lt2 € [d(l - 2e),d(l - 3s)), 5 , 3 G [d(l - 3s),d(l - As)), 5 1>3 G 
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[d(l — 4e),d(l — 5e)), and #2,3 S [d(l — 5e),d(l — 6e)). For each e there are £l(d 6 ) non-isomorphic 
matrices. If e is suitable small then all these matrices fulfill the triangle conditions and inequality 

LU □ 

In general we have a<(m,d) € ft(d m(m+1)/2 ), a<(m,d) € 0(d m ( m+1 ^ 2 ), a(m,d) € fi(d m ( m + 1 )/ 2 - 1 ), 
and a(m,d) e 0(d m ^ m+1 ^ 2 - 1 ). 

In Section [4] and Section [5] we give an algorithm to obtain a(d, 3) via implicit computer enu- 
meration. Some of these computed values are given in Table [l] and Table [2] For a complete list of 
a(d, 3) for d < 1000 we refer to [8]. This amounts to 

a<(1000, 3) = 4299974867606266 » 4.3 • 10 15 . 
3. Bounds for V 3 

In this section we give bounds for the probability V3 that A3 is the distance matrix of a tetra- 
hedron in the 3-dimensional Euclidean space E 3 , where A 3 is a symmetric 4x4 matrix with 
zeros on the main diagonal and the remaining entries being equi-distributed in [0,1]. Therefore 
we consider a discretization. Let d be a fixed number. We consider the d e six-dimensional cubes 
Cii,...,i 6 '■— x j=i [tp ^in — Pi I] 6 - For every cube C it is easy to decide whether every point of 
C fulfills the triangle conditions, no points of C fulfill the triangle conditions, or both cases occur. 
For inequality [l] we have no explicit test but we are able to compute a lower bound CMD^C) and 
an upper bound C M D^(C) , so that we have 



CMD 3 (C) < CMD 3 (x) < CMD 3 (C) for all x € C. 

Thus for some cubes C we can decide that all x E C correspond to a tetrahedron. We denote this 
case by S(C) = 1. If no x e C corresponds to a tetrahedron we set 5(C) = —1. In all other cases 
we define 5(C) = 0. With this we obtain for all d e N the following bounds: 

Lemma 4. 

C:H(C)=1 C:H(C)=-1 

Thus we have a method to obtain bounds on P3 using computer calculations. For the actual 
computation we use two further speed ups. We can take advantage of symmetries and use an 
adaptive strategy: We start with a small value of d and subdivide cubes C with 5(C) = recursively 
into 8 smaller cubes. After a computer calculation we obtain 

0.090 < V 3 < 0.111, 

which proves Theorem [2] Clearly Theorem [2] can be improved by simply letting the computers 
work for a longer time or by using a computing cluster, but the convergence of our approach seems 
to be rather slow. An enhanced check whether a cube C fulfills inequality ([!]) would be very useful. 

Good estimates for V 3 can be obtained by considering the values a(3,d) in the following way. 
At first we consider the probability V 3 being defined as V 3 where additionally <5o,i = 1. 

Lemma 5. 

V 3 = V 3 . 

Proof. The problem of determining Vz or V3 is an integration problem. Due to symmetry we only 
need to consider the domain where 5q,i is the maximum. For every <5o,i £ (0, 1] there is a probability 
p(fio,i) that 8o,x, ■ . ■ , 62,3 are distances of a tetrahedron where 60,2, ■ ■ • , ^2,3 €E (0, 80,1] are equi- 
distributed random variables. Since p(8q\) is constant we can conclude the stated equation. □ 
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Lemma 6. 

4-g(rf,3) 
V z = i im ' 

Proof. We consider a modified version of the algorithm described above to obtain exact bounds 
on 7^3 . As already mentioned, the triangle inequalities alone define a five-dimensional polyhedron. 
Since determinants are continuous CMD3 = defines a smooth surface and so the volume of all 
cubes C with 5(C) = converges to zero. Thus substituting 2(C) by the evaluation of S in an 
arbitrary corner of C yields the correct value for V3 — V3 for d —* 00. Since there are at most 
0(<i 4 ) six-tuples (d, i%, . . . , ig), ij € {1, . . . , d} with non-trivial automorphism group we can factor 
out symmetry and conclude the stated result. □ 

Using Lemma [2] and Theorem [I] we can modify this to 

_ aid. 3) ,, aid. 3) 17 

7> 3 = lim ' < lim ) ' = = 0.1416. 

d^oo a(d, 3) ~ d^oo a(d, 3) 120 

Heuristically we observe that the finite sequence ( — /j'ol ] is strictly decreasing. So the 

\ a l a ><V/ l<rf<1000 

following values might be seen as a good estimate for V%: 



a(600,3) _ 2000468396580 
a(600,3) ~ 19359502966749 * 

a(800,3) _ 8430487428682 
a(800,3) ~ 81665192828999 * 

a(1000,3) 25728695195597 



a(1000,3) 249377330461249 



0.103333, 
0.103232, and 
w 0.103172. 



4. Orderly generation of integral tetrahedra 

Our strategy to enumerate integral tetrahedra is to merge two triangles along a common side. 
In Figure [l] we have depicted the two possibilities in the plane to join two triangles (0, 1,2) and 
(0,1,3) along the side 01. If we rotate the triangle (0,1,3) in the 3-dimensional space from the 
position on the left in Figure [T] to the position on the right we obtain tetrahedra and the distance 
^2,3 forms an interval [l,u]. The restriction to integral tetrahedra is fairly easy. 




Figure 1. Joining two triangles. 



Let us consider the example 5q,i = 6, #0,2 = #1,2 = 5, £0,3 = 4, and 5i ( 3 = 3. Solv- 
ing CMD 3 = over the positive real numbers yields that the configuration is a tetrahedron iff 
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62,3 G ^ V702-24V45g ; ^702+^/455 j _ (2.297719304,5.806934304). Thus there are integral tetra- 
hedra for <5 2 ,3 € {3, 4, 5}. In general we denote such a set of tetrahedra by 

$0,1, $0,2, $1,2, $0,3, $1,3, $2,3 € [l,r]. 

This notation permits to implicitly list f2(d 6 ) integral tetrahedra in 0(d 5 ) time. 

All integral tetrahedra can be obtained in this manner. So an enumeration method is to loop 
over all suitable pairs of integral triangles and to combine them. We will go into detail in a while. 
Before that we have to face the fact that our enumeration method may construct pairs of isomorphic 
tetrahedra. Looking at Table [T] we see that storing all along the way constructed non-isomorphic 
integral tetrahedra in a hash table is infeasible. Here we use the concept of orderly generation [TS] 
which allows us to decide independently for each single constructed discrete structure if we have to 
take or to reject it. Therefore we have to define a canonical form of an integral tetrahedron. Here 
we say that a tetrahedron T with side lengths Sij is canonical if for the lexicographic ordering of 
vectors y, 

(#0,1, ^0,2, #1,2, $0,3, $1,3, $2,3) h ($r(0),r(l), ■ • • j K(2),t(3)) 

holds for all permutation r of the points 0,1,2,3. We describe the algorithmic treatment of a 
canonicity function x(T) 1— > {true, false} which decides whether a given integral tetrahedron T is 
canonical in Section [5] We have the following obvious lemma: 

Lemma 7. If x{$o,\, $0,2, 81,2, #0,3, $1,3, $2,3) = true and x{$o,x, £0,2, £1,2, £0,3, $\,3, $2,3 + 1) = false 
then x{$o,x, £0,2, £1,2, £0,3, £1,3, c>2 : 3 + n) = false for all n>l. 

Thus for given 6 ,i, 60,2, $1,2, 60,3, and £1,3 the possible values for <5 2 ,3 which correspond to a 
canonical tetrahedron form an interval [I, u). Clearly, the value of x(£o,i> £0,2, £1,2, $o,3, $1,3, $2,3) has 
to be evaluated for 5 2 , 3 G {5 itj -l,5i tj ,S itj + l | G {(0, 1), (0, 2), (1, 2), (0, 3), (1, 3)}} only Thus 
we can determine the interval [l,u] using O(l) evaluations of x(1~). 

Algorithm 1. Orderly generation of integral tetrahedra 
Input: Diameter d 

Output: A complete list of canonical integral tetrahedra with diameter d 
begin 

£0,1 = d 

for do, 2 from \_^-\ to d do 
for 5±,2 from d+ 1 — 60,2 to 60.2 do 
for <5o.3 from d + 1 — £0,2 to <5q,2 do 
for (5 13 from d + 1 — 60,3 to <5 .2 do 
Determine the interval [l,u] for 62,3 

print 60,1,60,2, £1,2, £0,3, £1,3, [I, u] 
end 
end 
end 
end 
end 

We leave it as an exercise for the reader to prove the correctness of Algorithm [T] (see [T5] for 
the necessary and sufficient conditions of an orderly generation algorithm]. Since we will see in the 
next section that we can perform x{T) in O(l) the runtime of Algorithm IT] is 0(d 4 ). By an obvious 
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modification Algorithm [T] returns a complete list of all canonical integral tetrahedra with diameter 
at most d in 0(d 5 ) time. 

We remark that we have implemented Algorithm [l] using Algorithm [2] for the canonicity check. 
For the computation of a(800, 3) our computer needs only 3.3 hours which is really fast compared 
to the nearly 3 hours needed for a simple loop from 1 to a(800,3) on the same machine. Due to 
the complexity of 0(d 4 ) for suitable large d the determination of a(d, 3) will need less time than 
the simple loop from 1 to a(d, 3). 

5. Canonicity check 

In the previous section we have used the canonicity check x(T) as a black box. The straight 
forward approach to implement it as an algorithm is to run over all permutations r G S4 and to 
check whether (Sq,i, ■■■ , 62,3) >z (#t(o),t(i)j • • ■ j ^t(2),t(3))- This clearly leads to running time O(l) 
but is too slow for our purpose. It may be implemented using 24 • 6 = 144 integer comparisons. 
Here we can use the fact that the tetrahedra are generated by Algorithm [T] So if we arrange the 
comparisons as in Algorithm [2] at most 6 integer comparisons suffice. 

Algorithm 2. Canonicity check for integral tetrahedra generated by Algorithm [T] 
Input: <5 ,i, £0,2, £1,2, ^0,3, £1,3, <$2 : 3 
Output: x(8o,u 80,2, £1,2, £0,3, £1,3, £2,3) 
begin 
if <5o,i = £0,2 then 
if £0,2 = 81.2 then 
if 5o,3 < £1,3 then return false 
else 

if (5i 3 < 62.3 then return false else return true 
end 
else 

if 5i t 3 < 62,3 then return false 
else 

if <5 ,i = #0,3 then 

if <5i 2 < £1,3 then return false else return true end 
else 

if <5o.i > <5i.3 then return true 
else 

if 6i t 2 < £0,3 then return false else return true end 
end 
end 
end 
end 
else 

if <5o.2 = £1,2 then 

if <5o.i < £2,3 or <5o,3 < <5i.3 then return false else return true end 
else 

if 80,2 = £1,3 then 

if £o,3 > £1.2 or £0,1 < ^2.3 then return false else return true end 
else 
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if #0,2 = #0,3 then 

if 5i,2 < #1,3 or <5o,i < #2,3 then return false else return irue end 
else 

if #0,3 > #1,2 then 

if #o,i < #2,3 then return true else return false end 
else 

if #o,i < #2,3 then return true else return false end 
end 
end 
end 
end 
end 
end 

To prove the correctness of Algorithm [2] we consider all vectors (#o,i, • • • , #2,3) with #jj € 
{0,...,5}, #0,1 > #0,2 > #1,2, #0,2 > #0,3, and #0,2 > #1,3- It suffices to show that Algorithm [2] 
returns the correct value for this finite set of vectors since these inequalities are fulfilled by Algo- 
rithm [l] and also necessary for x(T) = true. Algorithm [2] can be considered as a binary decision 
tree. It might be a task to optimize this type of binary decision tree in the worst or in the average 
case. 

6. Dimensions m > 4 

Clearly the question for bounds for P m arises also for m > 4. But non-trivial answers seem out 
of reach by our approach. So far we have no efficient equivalent of Algorithm [2] at hand and the 
number a(4, d) of integral 4-dimensional simplices with diameter d is fl (d 9 ) . We give the known 
values of a (4, d) in Table [6] 



d 


a(d,4) 


d 


a(d,4) 


d 


a(d,4) 


d 


a(d,4) 


1 


1 


14 


12957976 


27 


4716186332 


40 


162007000505 


2 


6 


15 


24015317 


28 


6541418450 


41 


202323976907 


3 


56 


16 


42810244 


29 


8970194384 


42 


251321436143 


4 


336 


17 


73793984 


30 


12168243592 


43 


310607982160 


5 


1840 


18 


123240964 


31 


16344856064 


44 


382002253424 


6 


7925 


19 


200260099 


32 


21748894367 


45 


467627887530 


7 


29183 


20 


317487746 


33 


28688094208 


46 


569910996879 


8 


91621 


21 


492199068 


34 


37529184064 


47 


691631229557 


9 


256546 


22 


747720800 


35 


48713293955 


48 


835911697430 


10 


648697 


23 


1115115145 


36 


62769489452 


49 


1006370948735 


11 


1508107 


24 


1634875673 


37 


80321260053 


50 


1207047969441 


12 


3267671 


25 


2360312092 


38 


102108730634 


51 


1442539675756 


13 


6679409 


26 


3358519981 


39 


128999562925 


52 


1718015775541 



Table 3. Number a(d, 4) of integral 4-dimcnsional simplices with diameter 1 < 



d < 52. 
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